Periodic Orbit Theory for Rydberg Atoms in External Fields 
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Although hydrogen in external fields is a paradigm for the application of periodic orbits and the 
Gutzwiller trace formula to a real system, the trace formula has never been applied successfully to 
other atomic species. We show that spectral fluctuations of general Rydberg atoms are given with 
remarkable precision by the addition of diffractive terms. Previously unknown features in atomic 
spectra are exposed: there are new modulations that are neither periodic orbits nor combinations 
of periodic orbits; 'core-shadowing' generally decreases primitive periodic orbit amplitudes but can 
also lead to increases. 
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Periodic orbit theory, in the form of the Gutzwiller 
trace formula (GTF) [Q, provides the most powerful 
framework for the semiclassical quantization of chaotic 
systems. It is more than a decade since it was first shown 
that the GTF provides a quantitative description of the 
oscillations in the density of states of highly excited hy- 
drogen atoms in magnetic fields j|. However, the trace 
formula has never been applied successfully to any other 
species of singly excited (Rydberg) atom. This is clear 
from comparisons between accurate quantal spectra since 
the spectral amplitudes for non-hydrogenic atoms differ 
substantially from those of hydrogen j| . 

Much effort has been expended in developing closed or- 
bit theory , the semiclassical theory which describes 
photoabsorption by atoms in external static fields. By 
matching semiclassical waves to Coulomb waves near the 
nucleus, the photoabsorption strength is obtained as a 
sum of contributions from only those orbits that close at 
the nucleus. In contrast, all periodic orbits contribute 
to the GTF (the eigenvalue spectrum). It was shown 
that, provided core-scattered waves are included consis- 
tently in the matching procedure, closed orbit theory 
can be applied to general atoms in fields But, 
for non-hydrogenic atoms one finds additional modu- 
lations of O(Vh) or higher relative to the hydrogenic 
contributions — the 'combination-recurrences' — that are 
due to sums of closed orbits. Contributions from the har- 
monics of closed orbits are reduced in amplitude through 
'core- shadowing' [§]-[l0| but those associated with the first 
traversals of the primitive orbits are unaffected. 

However, closed orbit theory does not account, even 
qualitatively, for the observed differences between atomic 
species for the eigenvalue spectrum, i.e. the density of 
states, p(E) = — Im Tr G(E) /it. In the density of states — 
quantitatively well-described by periodic orbit theory for 
hydrogen — amplitudes of primitive periodic orbits can 
vary substantially between atoms, in contrast with closed 
orbit theory. Also, the modulations associated with com- 
binations of periodic orbits appear at different orders in 
Ti relative to closed orbit theory. 

Rydberg atoms and molecules in the field-free case are 
described by Quantum Defect Theory (QDT), one of the 
most widely used theories in atomic physics. In QDT, 
the effects of a multi-electron core are described by a 



set of phase-shifts, or 'quantum defects', Si, in each par- 
tial wave, I. In the limit when the quantum defects vary 
smoothly with, I, they can be related to the classical pre- 
cession angle of the Kepler ellipse by 9 = 2irdSi/dl. For 
many atoms only the lowest partial waves have non-zero 
quantum defects. For example, for even parity lithium, 
Sq ~ OAn and Si>2 — while for helium Sq ~ 0.37T and 
5i>2 — 0. In this case, variation of Si with I is clearly not 
smooth. 

In this letter, we present an approach that, for the 
first time, combines the Gutzwiller trace formula with 
Quantum Defect Theory and hence sheds new insight 
on the classical interpretation of quantum defects. Our 
approach yields simple analytical expressions for the dif- 
ferences in amplitudes for general non-hydrogenic atoms. 
We compare the new semiclassical results with full quan- 
tal calculations and, for well isolated orbits, find them to 
be extremely accurate, for example to within about 1% 
at h ~ 0.01. 

A surprising finding, predicted by theory and con- 
firmed by quantal results, is that, although the am- 
plitudes of the primitive orbits are mostly reduced, as 
one would expect from the idea of 'core-shadowing', in 
the non-hydrogenic case they can also increase. This 
is shown below to be due to a de-phasing between 
diffractive (core) and geometric (Coulomb) contribu- 
tions. Combinations of periodic orbits appear with order 
at least h. Most significantly, new modulations appear 
that are not combinations of real periodic orbits but are 
rather pure diffractive orbits. They pass through the 
core and are made periodic by the diffraction. We em- 
phasize that all these effects are accurately described by 
the diffractive periodic orbit theory. 

The periodic orbit theory of diffraction was developed 
recently for Hamiltonians with discontinuities [|ll]-[l3). 
For our purposes, a good example of a diffractive sys- 
tem is the cardioid billiard, which has a single sharp ver- 
tex. In this case, periodic orbits are decomposed into two 
kinds: those that do not intersect the vertex (geometric 
orbits) and those that do (diffractive orbits). The density 
of states has been shown to be obtained as a sum: 



p(E) = ImTr G g (E) - -ImTr G D {E). (1) 

The first (geometric) term yields the ordinary GTF. The 
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trace, over the second (diffractive) contribution has been 
shown to be [OOl 



TyG d (E) = ?f II d(n)G(q n ,q n+1 ;E), 



(2) 



where T p is the total sum of periods taken over the paths 
between the vertices and d(n) is the diffraction constant 
which depends on the type of diffraction. Equation (|^) 
encapsulates the important result that the trace integral 
taken between the nth and n+lth vertices is proportional 
to the Green's function between those points. 

We apply diffractive periodic orbit theory to our 
atomic systems by treating the non-hydrogenic core as 
a diffractive source. The crucial step is to obtain an 
expression for the diffractive constant, d in terms of 
quantum defects. To this end, we consider an incoming 

Coulomb wave, V'coui' wmcn approaches the atomic core 
from infinity at an angle, Of, to the z-axis. On reaching 

the core, V'coui produces a scattered wave, V'scattj which 
feeds outgoing semiclassical waves along periodic orbits; 
V'scatt can be decomposed into an outgoing Coulomb wave 
together with a core-scattered wave || : ^ sca tt (r, 0) = 

^CouiM) + V'core (r,0). The Coulomb scattered wave is 
strongly back-focussed along ~ Of and can be written 
in closed form Q|. Our first approximation consists of 
equating V'ctmi with the source for geometric paths (i.e. 
the usual GTF). The core-scattered wave i/'core, arising 
from the incoming wave at angle Of, is equated with the 
source of diffractive semiclassical waves. At a radius, 
ro ~ 50 Bohr, we express "0core in a partial-wave expan- 
sion which, for m = is J^] 

Ac = 5- £ Y io(0f, 0)Yio(6, 0)e-^/ 4 (e 2 ^ - 1) 



1=0 



(3) 

where 5i are the quantum defects. Finally, we take d to 
be the fractional amplitude scattered by the core: 



d(0i, Of) = -0core(r O 



iVV'CoulO'O, 



(4) 



All calculations and comparisons with fully quantal 
spectra presented here have been carried out for s-wave 
scattering (appropriate for atoms such as lithium or he- 
lium which are used frequently in experiments of atoms 
in fields). So, below 5 = So and, in this case, ipcore is 
isotropic. However, generalization to odd parity or atoms 
with multiple quantum defects is straightforward. 

We consider the specific example of Rydberg atoms 
in a static magnetic field of strength, 7 (atomic units). 
The quantum spectra are calculated at constant scaled 
energy e = £7~ 2 / 3 , that is for fixed classical dynamics. 
Quantum mechanically, we calculate a set eigenvalues, 
7~ 2 ^ 3 , corresponding to different effective h [Q. Below, 
h denotes 7 1 / 3 . 

For Rydberg atoms in a magnetic field, the best stud- 
ied periodic orbits are the straight line orbit perpendic- 
ular to the field, and the 'balloon' orbit, V^. The 



well-known Garton-Tomkins orbit |lq] , Rx, is responsi- 
ble for the quasi-Landau oscillations observed in m = 1 
atomic spectra near the ionization limit at energy spac- 
ing ~ 1.5frui; these were the first observed 'footprints' of 
periodic orbits in a real physical system. The balloon 
orbit dominates m = 0, odd-l spectra with oscillations 
of spacing ~ 0.64?La;. The effect of the core on the orbit 
parallel to the field, V\, is relatively weak [p|p|9[]. The 
periodic orbit labelling terminology of Ref. |I6| is used 
throughout. 

For the case of s-wave scattering, each diffractive con- 
tribution in Eq. (0) is: 



dG = Vh(t 
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and, in effect, represents the contribution of a pure 
diffractive orbit. Note the additional phase of — ir/4 rel- 
ative to an equivalent geometric primitive periodic orbit. 

Now we see that the amplitude of each non-hydrogenic 
primitive periodic orbit actually arises from the inter- 
ference between two contributions with the same action 
but different phase: a geometric one of the Gutzwiller 
form weighted by the trace of the stability matrix, M, 
in the usual manner, i.e. for the pth orbit = 
(ttS p ) I W 1 2 — Tr M p \ , and a diffractive one, given by 
Eq. (^), following the same path and of similar action 
but weighted by lf-y/rrivi, where mi2 is an element of M. 
This contrasts with the cardioid billiard where a typical 
contribution is either pure geometric or pure diffractive. 

We can easily show that the fractional reduction of am- 
plitude for a primitive periodic orbit of a non-hydrogenic 
atom relative to that of hydrogen is 

A s /A H = \Jl + AK 2 sin 2 (5) h-4K sm(5) sin((5 + 0) Vh, 

(6) 

where, in general, 



S27T 

■R = A H \j sin(6» i /2)sin(6»//2). 

77742 



(7) 



All parameters, e.g. initial and final angles, 0i, Of, re- 
fer to the particular primitive orbit under consideration. 
The fractional reduction or increase is oiO(Vn). In gen- 
eral, <p = ~ 7r/4. However, R± runs along a boundary 
of the fundamental symmetry domain and so requires 
special treatment: when stable, with winding number v, 
we find <f> — — 7r/4 — vi: /2 and 1Z = 2^J(tx jm\2) sin(-7ri/). 
The most important correction in Eq. (^J) is the \ffi term. 
This is zero for 5 = —<j) and positive for sin(<5 + </>)< 0, 
leading to an increased amplitude for a primitive periodic 
orbit. In contrast, in closed orbit theory, the main pho- 
toabsorption source term and the core scattered terms 
do not have this 7r/4 de-phasing. 

Harmonics of primitive orbits also have further contri- 
butions from product terms of O(h). However, product 
terms also give additional weak contributions at actions 
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that are sums of periodic orbits. For example a two-orbit 
combination has amplitude 



8^ 2 (Si + S 2 ) 
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(8) 
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FIG. 1. (a) Comparison of Fourier transforms of the den- 
sity of states for hydrogen and lithium (S = 0.40 457r) in a 
static magnetic field at constant scaled energy e — —0.2 from 
a fully quantal calculation with average h = 1/90. Note the 
changes in amplitudes of periodic orbits and new modulations 
due to diffractive orbits in the lithium case, (b) Comparison 
between quantal and semiclassical difference spectra obtained 
by coherently subtracting the Fourier transforms shown in 
(a). This exposes the diffractive contributions to the spec- 
trum and eliminates contributions from orbits which do not 
pass through the core. Shown are changes in periodic orbit 
amplitudes due to diffraction, diffractive combinations of two 
periodic orbits and pure diffractive orbits, marked Di and 
I>2. Away from bifurcations, which affect Vi and D2, the 
agreement between quantum and semiclassical calculations is 
excellent. 

In Fig. [l](a) we show Fourier transforms of the oscil- 
latory part of the even-^, m = 0, eigenvalue spectra for 
hydrogen and lithium (S = 0.40457r) in a static mag- 
netic field at constant scaled energy e = —0.2 and with 
n = 7- 1 / 3 = TiT 1 ranging from 60 to 120. In Fig. 0(b) 
we plot the 'difference' spectrum obtained by coherently 
subtracting the Fourier transform of the hydrogenic spec- 
trum from that of lithium; this exposes the diffractive 
contributions and eliminates contributions from periodic 
orbits which do not pass through the core. For com- 
parison, we also plot a semiclassical difference spectrum 
obtained by summing all terms of order y/h and h; agree- 
ment is excellent. The discrepancy in D2 is due to the 
effects of bifurcations that are not taken into account in 



the semiclassical calculation presented here. 

We can see that for lithium the amplitudes of Ri and 
its harmonic R2, as well as V] 1 and other orbits are 
substantially reduced. There are additional small peaks 
which correspond accurately to sums of periodic orbits. 
Importantly, there are strong peaks (marked D\ and D2) 
which do not match any combination of orbits. At these 
scaled actions (S ~ 2.87 and S ~ 2.94) we find orbits that 
are closed but not periodic. For hydrogen, only orbits 
that are periodic in the fundamental symmetry domain 
contribute. Here we see that pure diffractive orbits, such 
as D 2 , can contribute to the non- hydrogenic spectrum at 
0(y/K) so are substantially stronger than combinations 
of orbits. The peak at S ~ 2.87 is due to an isolated 
closed orbit and is obtained almost exactly from Eq. (||) 
as seen in Fig. |l|(b) (note that in Fig. |l|(a) the peak asso- 
ciated with this orbit is masked by the peak of a periodic 
orbit which does not approach the nucleus). The peak 
at S ~ 2.94 consists of contributions from a pair of non- 
isolated orbits close to a bifurcation so their contribution 
is over-estimated semiclassically. On examination of the 
diffractive orbits we find that they correspond to the first 
closure of asymmetric periodic orbits, some of which cor- 
respond to the X n series of 'exotic orbits' In hy- 
drogenic eigenvalue spectra such orbits only contribute 
at their full period, whereas in the diffractive case they 
appear at closure. 

We have carried out a detailed study of these effects for 
several scaled energies to study the Ti and 6 dependence 
of the diffractive effects. In Figs |^(a)-(d) we compare 
the fractional change relative to hydrogen between the 
fully quantal and semiclassical expressions for R\ and 
V/ . The agreement is very good. For the h dependence 
there are fewer points for V] 1 since a wide spectral range 
is required to resolve it from a nearby orbit. An especially 
interesting feature is the de-phasing of i?i relative to V\ 
seen in Figs ||(c) and ||(d) . The diffractive contribution 
to Vi is — 7r/4 out of phase with the geometric term. As a 
result the amplitude exceeds that of hydrogen for <5 < 7r/4 
and is minimal at S ~ 0.657T. In contrast, the geometric 
and diffractive contributions for i?i are almost in phase 
at e = —0.275 and remain so for a wide range of scaled 
energies about e ~ —0.3, where the orbit undergoes its 
2 : 1 resonance with v ~ 0.5. 

In Figs ||(e) and ^(f) we investigate the combination 
orbits and the diffractive orbit that appears at e = —0.2 
for S ~ 2.87. In this case we plot the ratio of amplitudes 
relative to the first peak of C, the circular orbit, a pe- 
riodic orbit that does not pass through the nucleus and 
hence is unaffected by the diffraction. In both cases the 
agreement is very good. 

In conclusion, we have shown that periodic orbit theory 
(the GTF) may be applied to all singly excited atoms as 
successfully as for hydrogen by bringing in the effects of 
QDT in the form of diffractive corrections. Also, to our 
knowledge, this is the first demonstration of a diffractive 
effect in a real system, since previously diffraction has 
only been applied to model problems such as billiards. 
Although we show explicit results for s-wave scattering 
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in lithium and helium, our method is applicable generally 
to other Rydberg atoms and molecules in external fields. 
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FIG. 2. Dependence of diffractive contributions on h and 
8: comparison between quantum results (full circles) and 
semiclassical formula, Eq. ([]), (solid line) showing near ex- 
act agreement. The vertical axis represents the ratio of 
non-hydrogenic to hydrogenic amplitudes, (a) Dependence 
of As/ A H on h for ft at e = -0.2, -0.6 and -0.45. (b) 
Same as (a) for V\ at e = —0.2, and —0.275. (c) Dependence 
of Ag I Ah on quantum defect, 8, for Ri. (d) Same as (c) for 
Vi : note the de-phasing relative to Ri and that the ampli- 
tude exceeds the hydrogenic value for 8 < 7r/4. (e) and (f) 
dependence of contributions of combination and diffractive 
orbits on 8. Here the ratio relative to the first contribution of 
the circular orbit, C, is shown. As predicted by the theory, 
combination orbits show a sin 2 5 behaviour while diffractive 
orbits follow a sin 8 curve. 

An interesting recent calculation JllJ treated the pho- 
toabsorption of general atoms within the framework of 
the standard theory using a model potential. Then, the 
observed closed orbit modulations were modelled by su- 
perposing thousands of very unstable orbits. Hence, the 
issue of whether the dynamics of non-hydrogenic atoms 
at moderate scaled energies is an instance of chaos (i.e. 
very unstable motion) or an effect 'beyond periodic or- 
bits', such as diffraction, remains open. Our work ad- 
dressed this issue. 

Currently, there is added interest in diffractive systems 
since they have very recently been associated with a new 
class of intermediate level statistics ('Half-Poisson') Jl8| . 
Recently, eigenvalue statistics for rubidium in fields 
were investigated experimentally |l9[ ] and shown to be 
nearer the GOE limit than comparable hydrogenic re- 
sults. Hence, given that spectroscopic resolution exceeds 
mean level spacing, experimental verification of diffrac- 
tive effects in the eigenvalue spectrum, for example the 
presence of the 'D' modulations, is in principle possible. 
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